Three-body properties in hot and dense nuclear matter 
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We derive three-body equations valid at finite densities and temperatures. These are based on 
the cluster mean field approach consistently including proper self energy corrections and the Pauli 
blocking. As an application we investigate the binding energies of triton and determine the Mott 
densities and momenta relevant for a many particle description of nuclear matter in a generalized 
Beth-Uhlenbeck approach. The method, however is not restricted to nuclear physics problems but 
may also be relevant, e.g., to treat three-particle correlations in weekly doped semiconducters or 
strongly coupled dense plasmas. 
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I. INTRODUCTION 

Correlated many particle systems such as, e.g., nu- 
clear matter or strongly coupled plasmas, have a compli- 
cated dynamical behavior. Only few areas of the density- 
temperature phase diagram used to characterize the state 
of the system in thermal equilibrium may be described in 
the approximation of noninteracting quasiparticles. The 
dynamics of the quasiparticle is determined by the mean 
field of the other particles and therefore some promi- 
nent features like self energy corrections and Pauli block- 
ing (Bose enhancement) are sufficiently taken into ac- 
count ||l|. However due to sizable residual interactions 
many interesting and exciting phenomena, such as clus- 
tering, formation of condensates, and phase transitions 
occur. 

A large number of these phenomena such as, e.g., two- 
particle bound state formation, can already be accounted 
for by explicitly introducing two-body correlations into 
the formalism. This may be achieved in the frame work 
of the Green function method M and leads to effective 
two-body equations that include medium effects in a con- 
sistent way. 

Recently, generic three-body processes have been cal- 
culated utilizing exact few-body methods in the context 
of many particle systems for i) nuclear matter P,p| and 
ii) plasmas at star conditions [0 . Although the relevance 
of three-body processes in many particle systems have al- 
ready been recognized, little progress has been achieved 
since both fields (many-particle and few-body physics) in 
itself are rather elaborated, bearing their own technical 
problems. However, some phenomena of many particle 
systems require a treatment of effective few-body sys- 
tems embedded in a medium, in particular three-body 
processes. 

In the context of nuclear matter, a three-body 
Faddeev-type equation has been derived within the Green 
functions method applicable to describe three-body cor- 



relations in nuclear matter of thermal equilibrium P,pj . 
Empirical evidence, including recent experimental data 
on cluster formation |H,pl, indicate that a large frac- 
tion of deuterons can be formed in heavy-ion collisions 
of energies below E/A < 200 MeV. The abundances of 
deuterons are determined by the deuteron formation via 
NNN -^ dN {N nuclcon, d deuteron) and break-up, 
dN -^ NNN, reactions. The hadronic reaction requires 
a proper treatment of the effective three-body problem. 
We have numerically solved the scattering problem to 
describe the deuteron break-up reaction consistently in- 
cluding the effects of finite densities and temperatures of 
the surrounding nuclear matter. Within linear response 
theory we have given a first estimate for the break-up 
time of deuterons embedded in nuclear matter g] . 

Indeed, at moderate densities n ^ f^o/10 (no = 
0.17 fm^ nuclear matter density) and temperatures T < 
15 MeV nuclear matter may be considered as a mixture 
of light nuclei (chemical picture). The transition into 
that region is determined by the Mott density M and 
is relevant, e.g. for the final stage of a heavy ion colli- 
sion at intermediate energies, respectively on the surface 
of the expanding nuclear matter. The Mott effect has 
also been considered at the beginning of the heavy ion 
reaction using inverse photo disintegration pi. Here we 
present a calculation that determine the Mott densities 
for the three-body nuclear bound state. 



II. BETH-UHLENBECK APPROACH TO 
NUCLEAR MATTER 

The basis to treat correlated densities is provided by a 
generalization of the Beth-Uhlenbeck approach ||l3] . The 
nuclear density n — n{^,T) as a function of the chemi- 
cal potential fj, (for the time being we assume symmetric 
nuclear matter) and temperature T may be written as 



n = nfrec + n-c 



rih. 



^Y.h 



(1) 



For abbreviation we use 1 = ai = {fci, si,ti, . . .} de- 
noting momenta, spin, isospin, etc. of particle 1. The 
one-particle Fermi function is 



/i = /(ei) = (exp[/3(ei-M)] + !)-!, 



(2) 



where ei, denotes the quasiparticle energy and (3 the in- 
verse temperature. In first iteration the correlations may 
be treated on the basis of residual interactions between 
the quasiparticles. To do so the imaginary part of the 
self energy I](l,a;) should be small and the respective 
spectral function may be expanded with respect to the 
imaginary part S/(l,cj) as explained in Ref. pOJ] . The 
contribution of the correlated density to the total den- 
sity may then be written as pO[ 

ncorrifJ., ^) = X! / 27r ^^'-^' ^ ^ ^'^^^ 



X [/H-/(£i)] 



V 



duo Si — cj ' 



(3) 



where V denotes principle part integration. To evaluate 
the self energy we utilize cluster decomposition taking 
into account two- and three-particle, in general n-particle 
correlations (see Fig. |l|), 

E(l, z,) = Y,Y. ^2(12, 12; z, -f zx)Gi (2, za) 

2 A 

+ Y.Y. ^3 (123, 123; z, + ZA + zy) 



2,3 AA' 



xGi(2,za)Gi(3,za') 



(4) 



where Gi{zy) denotes the one-particle Green function of 
Matsubara frequency z^ H] , and the index c denotes con- 
nected terms only. The n-body t-matrices are defined in 
terms of the n-body Green functions G„(z) via 



G„(2) = Gl")(z) + Gi°)(z)r„(z)Gl")(z) 



(5) 



Using the spectral representation of the t-matrices and 
evaluating the Matsubara sums eventually leads to the 
following decomposition of the correlated density induced 
by the cluster decompostion Eq. (0), 



2n2 -I- 3n3 



"-2 



'2 1 ■ 



(6) 



where 712(3) denotes the two- (three-) particle correlated 
densities, present as bound n° or scattering rf^ states 
in chemical equilibrium. For the two-particle case the 
correlated densities have been given explicitly in Ref. fin] 



P>PMatt 



3 g{Econt + Eh), 



E ^9iEc. 

P>PMatt 



E 



duj 

2^ 



dco 



g{E cont + L^) 



y Ca ■ {Sa — sin Sa COS 6a ) , 



(7) 



where the sum over a indicates partial wave decompo- 
sition with Ca the proper Clebsch-Gordan coefficient and 
Sa the corresponding phase shift. The energies appearing 
in Eq. (R) are the binding energy i?b and the continuum 
energy i?cont that includes the proper self energy shifts. 
The Bose function for the two-fermion system is 



g{uj) = (e 



- Cp/3{"-2m) _ n-i 



1)- 



(8) 



To evaluate the contribution of the three-body correlated 
density n^ we presently focus on the bound state contri- 
bution. The i-matrix is then dominated by the bound 
state pole of energy E\^ . Evaluating the respective Mat- 
subara sums for the three-body bound state contribution 
in Eq. (^ using Eq. (Q)and the spectral representation of 
T3 leads to a rather simple expression for the three-body 
bound state contribution rtg to the total density, viz. 



E 

-P>PMott 



4(e 



/3(Boo„t+Bb-3M) 



(9) 



where i?cont the respective continuum energy. 

At moderate densities not only two-particle but also 
three- and four-particle correlations occur. Here we con- 
sider the three-nucleon bound state within the general- 
ized Faddeev approach discussed in the next section. One 
obvious feature in the above equation is the appearing of 
the Mott momentum PMott • It is important in the frame- 
work of the cluster Hartree-Fock expansion MKU closely 
related to the self consistent RPA [|l2[ and extended to 
finite temperature in |3,H.n3 that are followed here to 
consistently treat few-body correlations in matter. 



III. FINITE TEMPERATURE THREE-BODY 
EQUATIONS - BOUND STATES 

The formalism to derive few-body Green functions 
within the cluster mean field approximation at finite tem- 
peratures and densities has been given elsewhere P,[13[. 
Here we give some of the basic results and extend the 
formalism to include bound states. The reaction cross 
section based on the suitably modified AGS formalism 
has been given in [|| and the life time of deuteron fluc- 
tuations in nuclear matter as a first application in 0. 

The hierarchy of Green function equations is truncated 
using cluster mean field expansion, and a calculable form 
is achieved by introducing ladder approximation. We 
consider generic elementary two-particle interactions V2 
only. The resulting equations for the decoupled one-, 



two- and three-body Green functions at finite tempera- 
tures (utilizing the Matsubara technique to treat finite 
temperatures) will be given in the following. The one- 
particle Green function reads 



Gi{z) = r["\z) = {z - ei)-\ 



(10) 



In mean field approximation the single quasiparticle en- 
ergy £i is given by 

Zmi 

^^""(1) = E [^2(12, 12) - ^2(12, 21)] h, (11) 



and /2 — f{^2) the Fermi function given in the previous 
section. The equation for the two-particle Green function 
G2{z) reads 

G2{z) = N2R'i\z) -f Rf\z)N2V2G2{z), (12) 

where the two-body resolvent i?2 {z) is given by 



811,82 



4°^(12,l'2';z) 

Z - £1 - £2 

and the Pauli blocking factor A'2 by 



(13) 



N2{12, 1'2') = <5n„522'533'(/i/2 - /1/2) 

= 5ll'(522'fe'(l-/l-/2)- (14) 

We use the notation / = 1 — /• The respective equation 
for the three-particle Green function relevant to describe 
three-body correlations in a medium is given by 

G^{z) = N^Rf{z) + Rf{z)W^G^{z), (15) 

where the effective potential W3 reads 
3 



W3{l2?,,l'2'-i') ^^W^i\l2?,,l'2'i') 



(16) 



fe=i 



^(3)/-ioQ 1'0'q' 



W^'> (123, 1'2'3') = (1 - /i - /2) V^2 (12, 1'2')<533' • (17) 

The last Eq. is given for fc = 3 and cyclic permutation is 
understood. Note that W3. ^ VFJ. The Pauh factors Nz 
and the resolvents -R3 read respectively 

7V3(123, 1'2'3') = <5ii'<522'<533'(/i/2/3 + /1/2/3) (18) 

i?f (123,l'2'3';z) == ^""^^^'^^^' . (19) 

■^ ' z - £1 - £2 - £3 

Note that [A'3, R\^'] — 0. For convenience we may intro- 
duce the Green function of the nonintcracting system, 
viz. 



.(0) 



Gt'\z)^N,Rt>{z). 



(0), 



(20) 



If we now introduce a potential V3 — N^ W3 we may 
instead of Eq. ([l5|) write 



.(0) 



G3{z)^G'i"iz) + G'i'>iz)V3G3{z), 



(0), 



(21) 



which looks formally as the equation for the isolated 
case [|lj| and allows one to use three-body techniques to 
arrive at a solvable form, e.g. AGS equations for the tran- 
sition operator [pHl5||. These will be used to derive nu- 
merically solvable Faddeev type equations at finite tem- 
perature. Although we assume a Fermionic system the 
proper symmetrization is treated separately. 

Already in Eq. (|6|) we have introduced the channel no- 
tation that is convenient to treat systems with more than 
two particles [|4|. In the three-particle system usually 
the index of the spectator particle is used to characterize 
the channel. 

If the correlated pair, e.g. (12) and the spectator par- 
ticle, e.g. 3, are uncorrelated in the channel (3) we may 
define a channel Green function G3 (z). We generalize 
to the channel (7). The channel Green function is then 
defined by 

^3''W = zb^ ^G2(^a)Gi(z-^a). (22) 

The summation is done over the Bosonic Matsubara fre- 
quencies uj\, A even, lox — 'KX/{—il3) + 2/i. The equation 
for the channel Green function is derived in the same 
way as for the total three-particle Green function given 
in Eqs. ^ and (|l|). The result is 



d^\z) = Gf\z) + Gf\z)vt^Gi-^\z), (23) 



(7) _ 



(no summation of 7). Introducing the notation K 

V3 — V3 we arrive at the following equation for G3(z) 

expressed through the channel Green functions G3 (z). 



-.(7)/ 



G3(z) = G^-^^(z) + G^-^^(z)1^3^-^G3(z). 



1(7)/ 



.(7), 



(24) 



Now we have set the necessary equations, i.e. 
Eqs. (21), ( p3| ) and (EJ) to derive a proper integral equa- 
tion for scattering (see Ref. [^i||) and bound states at 
finite temperatures and densities. 

To arrive at a homogeneous three-body equation for 
the bound state |^() in medium we insert Eq. ( pi|) into 
the Lippmann-Schwinger equation 



|^t) = lim zeG3(^t+ze)|*t) 

£—►0 

After performing the limit we find 

\^t)^Gf{Et)Vz\^t). 
The Faddeev components are given by 

|^(")) = Gf ^3^"^ \^t) 



(25) 



(26) 



(27) 



and 



i*t) = Ei*^"^ 



(28) 



To arrive at a conveniently solvable version of the three- 
body bound state problem in medium we introduce form 
factors. To this end Eq. (p4) is inserted into Eq. (25). 
The limes e — > results in 



|*,)=G("'(i?OV^i"'|*t). 
Introducing the usual form factors 

eventually leads to an integral equation 

P 



(29) 
(30) 
(31) 



and finally to the bound state given in terms of the form 
factors 



|.p,)..^GfTi«Gf |f(« 



(/9) 



The transition channel operator Tg is defined via 



Gf' 



Gf 



^(0)^(/5)^(0) 



(32) 



(33) 



inserting this equation into Eq. (|2^) leads to an equation 
for the channel t matrix 



T, 



(/3) 



.(/3) 



K'^^ + crK'^r, 



iO).AP)rj.i0) 



3 ' 



and to V^^^Gi^^ 



T^'^^G^°\ Writing the PauH factors 
explicitly in the integral equation ( pl| ) results in 

|^W)^^(l_<5„^)rf)iV3i?(o)|F('^)). (35) 

P 



Note that G 



(0) 



duce t: 



■(/3) 



A^3i?f 



i?3 N3. We may now intro- 



N^/^T^^^N^/^ and |f*(«) = N^/^\F(P^). 
The resulting equation is 

|i.*(")) ^^(l_j^^)T3*('3)i?f |^*(«). (36) 

The equation for the transition channel operator T3 is 

then(V^3*(^)=iV3^/V3"''7V3^/^) 



where we have used (/1/2/3 + /1/2/3) = (1 - /» - /j)(l - 
fs + 9(^1 + £2)) valid for all permutations of ijk = 123. 
Note that 1 — /i — /2 > which is the case for low den- 
sities and the last equality in Eq. ( |38| ) holds for /^ <S /. 
Utilizing this approximation the corresponding scattering 
solution using a separable ansatz for the strong nucleon- 
nucleon potential has been given in H] . In Ref. H we have 
calculated the break up cross section Nd -^ NNN and 
found considerable dependence of the deuteron fluctua- 
tion time on the proper treatment of the medium depen- 
dence 1^. Here we solve Eq. ( p^ ) for the Yamaguchi p7| ] 
and a high rank separable version (Paris (EST)) of the 
Paris potential jlSJ . For a detailed overview on the pro- 
cedure we refer to [|l^. To compared to the perturba- 
tion theory result we use the respective wave function of 
the isolated triton and evaluate the medium dependent 
part, i.e. the Pauh blocking — (/i + /2)V2 and self energy 
corrections Y,^^{1) of the effective Hamiltonian in the 
standard fashion. 



IV. KINEMATICS 

Unlike the isolated three-body problem Galilei invari- 
ance (for the three-body system embedded in a medium) 
is not satisfied, since the Fermi functions depend explic- 
itly on the relative momentum of the three-nucleon sys- 
tem with respect to the medium. 

As a consequence one has to solve the three-body 
problem at a finite center of mass momentum Pcm = 
fci + ^2 + ^3 in a medium that may be considered at rest 
(34) ^mcd = for simplicity, see Fig. ||. However, techni- 
cally it is more convenient to let the three-nucleon sys- 
tem rest Pern — and the surrounding medium move 
with Pnicd = — (fci -I- ^2 -I- ks). This procedure results in 
the least change of the three-body algebra and is possible 
since the dependence on Pcm — Pmed is only through the 
Pauli blocking factors thus parametric. 

For simplicity we use angle averaged Fermi function 
((N2)) and ((7V3)), i.e. 



{{Ns)) 



1 



(47r)2 



d cos 61, d cos 9p d(j)qd4ip iV3(p, q, P^n), 



(39) 



where the angles are taken with respect to Pcm and p, q 
are the standard Jacobi coordinates iflJ. 



„*(/3) _ -,,*(/3) 
-^3 — ^3 



V, 



*(/3)^(0)y*(/3) 



(37) 



V. RESULTS 



Inserting all definitions the explicit form of the effective 
potential arising in this equation reads 

F3*(^)(123,l'2'3') = iV2'^'(12)(l - f3+9iei+e2))-'/' 
xF2(12, 1'2')<533'(1 - /3 + giev + S2')y/^N',/\r2') 



(1 - /i - /2)'/'V^2(12, 1'2')(1 - fv - /2')' 



/2 



(38) 



Since the Green functions have been evaluated in an 
independent particle basis the one-, two-, and three- 
particle Green functions are decoupled in hierarchy, as 
given in Eqs. ([l0|), (|l2|), (^. To solve the in medium 
problem up to three-particle clusters the one-, two- and 
three-particle problems are consistently solved. This 



leads to the single particle self energy shift Eq. (11), 
the two-body input including the proper Pauli blocking 
Eq. (1371) , and eventually to the three-body bound state 
Eqs. (|32|) and (|3q ) (or scattering state). 

For technical reasons we have ap pro ximated the nu- 
cleon self energy calculated via Eq. { \\/\\) by use of effec- 
tive masses for the nucleon, which is reasonable for the 
small nuclear densities considered here. 

For comparison we have calculated the energy shift of 
the three-body bound state using perturbation theory 
and the solutions of the respective isolated system. 

The resulting binding energies per nucleon in case 
where triton and medium are both at rest, i.e. Pcm — 
Pmcd = 0, as a function of the nuclear uncorrelated den- 
sity n are shown in Fig. ^ for two different temperatures 
T = 10, 20 MeV. The rank one Yamaguchi potential has 
been used for the solid lines, also used earlier in the con- 
text of in medium break up reactions {Nd — » NNN) [g . 
The dashed line originating at the deuteron binding en- 
ergy Ed/ 2 ~ —1.11 MeV reflects the respective contin- 
uum threshold (since the deuteron also changes the bind- 
ing energy the continuum threshold changes). The inter- 
section between the triton binding energy and the con- 
tinuum defines the Mott transition density. The Mott 
transition of the triton leads for both temperatures con- 
sidered directly to a three-body break-up. 

For T = 20 MeV the long dashed line shows the result 
of the Paris (EST) potential. The difference between the 
dashed line and the respective result of the simple Yam- 
aguchi potential is mostly due to the difference in bind- 
ing energy. The shape of the curves, also for the case 
T = 10 and 30 MeV (not shown) are very similar. The 
dashed dotted line is the result of a perturbative calcu- 
lation of all medium effects evaluated by using the triton 
wave functions for the isolated case. The dominant effect 
comes from Pauli blocking. The self energy correction in- 
cluding the continuum shift (due to effective masses) is 
smaller. 

The momentum dependence of the binding energy per 
nucleon for T = 10 MeV is shown in Fig. ^. As expected 
and known from the deuteron case the influence of the 
surrounding matter on the binding energy is decreasing 
with increasing momentum. Therefore the Mott transi- 
tion moves to higher densities for higher momenta. As 
a consequence in a moving system (like in a heavy ion 
collision) deuterons, tritons and presumably other light 
clusters may be formed at higher densities than expected 
from the simple considerations at rest. The momentum 
where the triton binding energy crosses the continuum is 
referred to as Mott momentum Piviott- 

Finally, Fig. |5| shows the dependence of the Mott mo- 
mentum on the density for the triton and the deuteron. 
The momenta and the densities related to the area above 
the respective curves allow bound states to be formed (i.e. 
Et,d < 0). At higher densities (respectively at higher mo- 
menta) the temperature dependence is less pronounced 
than for the lower densities. 



VI. SUMMARY AND CONCLUSIONS 

We have given proper equations to treat three-body 
correlations in matter and solved them for the nuclear 
bound state. The framework is provided by the cluster 
mean field (cluster Hartree-Fock) approximation closely 
related to the self consistent RPA approach, known for 
T = 0, and generahzed to T ^ recently |,^. 

We find substantial changes of the triton properties 
(being a much simpler problem the deuteron has been 
discussed much earlier [g|) and give the respective Mott 
densities and momenta. 

The application of rigorous few-body methods in the 
context of a many-particle description of nuclear matter 
provides a fruitful method not only useful for nuclear 
physics. Potential areas for a new application may be, 
e.g., weekly doped semiconductors to include the effects 
of bound eeh charged exciton states (trions) or strongly 
coupled dense plasmas where ionisation rates may rely 
on few-body reactions if the charge of the ions become 
large. 

The method is capable to be applied to the a par- 
ticle as well, which is the strongest bound nucleus and 
therefore very important for nuclear matter. Recently, 
indication for a-particle quartetting (condensation) has 
been found in a rather simple approach that need to be 
confirmed [Eol. 
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FIG. 1. Cluster decomposition of the self energy E in 
terms of n-body t-matrices. The sum is over irreducible con- 
tributions only. 




FIG. 2. Kinematical variables of the three-body system 
embedded in a medium. 




FIG. 3. Triton binding energy per nucleon (solid lines, 
Yamaguchi potential) as a function of the uncorrelated nu- 
clear density at a given temperature T; long-dashed lines: 
corresponding pertubation result, dashed-dotted lines: Paris 
(EST) potential. Dashed lines show the Nd continuum 
threshold for T = 10 MeV (left) and T = 20 MeV (right). 
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FIG. 4. Triton binding energy as a function of nuclear 
density at T = 10 MeV, and total momentum relative to the 
medium. From left to right: P = 0, 1, 2, 3, 4, 5 fm~^. 
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FIG. 5. Mott momentum per nucleon Pmou/A for triton 
- solid (T = 10 MeV) and dashed line (T = 20 MeV) and 
deuteron - dashed dotted (T = 10 MeV) and short dashed 
line (T — 20 MeV) as a function of nuclear density. Below 
the respective lines tritons (dcutorons) do not exist as bound 
states. 



